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La presente invention \a pour objet une niethpde pour f fprmer un mpdule a reseaux 
5 neuronaux, destine a^simuler en temps reel le mode d'ecbulemen~fr en tout 1 point d'une 
cpnduite, d'une veirie de fluidej>olyphasjq^ au moiris uhe/jphase liquide s et au 

/moms une phase gazeuse,' qui : soir le mieux adapte a des,. conditions op v eratoires et a viin 
/ • ensemble de grandeurs:physiques;'fl«ees -\ l r *y "V- t . \- \ A 

/• ■ - m 'tM<? \\\ v 

,Etat de la technique \'P:h M^'H^T'' v '\\ 
10 ^GSE' acheniinement des hydrdc$r$u^^ de production^ Qsqu'aifeutiites 




de traitenierit/constitue un mailldnlii^ortaht^cle latfelTaine;>petroliere. G'est un maillon 
^"delicat en raison de la complexite dfes'inte 4 rMt^6^'lferih*e " lS^s fpftases constituahtUes effluents 
transported >Les operateurS vOnt^pouriSbje^ d'|atieindre\une product! vate optimale 

S\ ill ^Jllfr 1 :^^: p) LU 

dans lessmfeilleures-'eoriditidnjs de securite^Bs^doiyent doriG*gerer au mieuxMajvitesse eUa 



/a/modelisfer'au 



, . y^) 

\ tmfeux le transpMiM!/-fi^ instiantjune 
\ image des ecoulerfiehts dans les^differentes pafties^de la chaine/de>productionixtenant 
* \. compte**de la constitution precise de l'effluent, les debits expressions etjTgs^ modes 
^eco^enj ^^^^^ ^ ^^0^^ ^ ^ 

^\\I1 exist^actuellement differents^m flux 
polyphasiques complexes r permettant a une stade precoce de . epnctevoirldes equipements de 



production adages. 



°fc LA 1 ^ 



Par les breve^U^gOTei, FR 2.756.044 (US6^^2')"et FR 2 756 045 (US 5 
25 960 187) du demandeur, notamn^t^rDconMti^des-'modules de modelisation formant le 
modele TACITE connu des gens de metier, permettant de simuler le transport de flux 
polyphasiques complexes en regime permanent ou transitoire et capables de prendre en 
compte des phenomenes d' instability qui se produisent du fait de la geometrie irreguliere 
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du terrain ou passe la conduite ou de la topographie de celle-ci, que les specialistes 
designent par « terrain slugging » ou « severe slugging ». 

La complexite des modules de simulation est a 1' image de celle des ph6nomenes 
modelises. Pr6cision et performances ne peuvent 6tre obtenues qu'apres un temps de 

5 moderation relativement important qui s'avere difficilement compatible avec une gestion 
en temps reel. C'est pourquoi ces modules de modelisation ne peuvent en l'&at servir a la 
gestion en temps reel de la production. Dans cette optique, il est ainsi apparu necessaire de 
mettre en ceuvre des methodes de modelisation permettant un bon compromis entre rapidite 
de calcul et precision des resultats. 

10 Par la demande de brevet FR 00/09.889 du demandeur, on connait une methode pour 

estimer en temps reel le mode d'ecoulement en tout point d'une conduite de structure que 
Ton peut definir par un certain nombre de parametres de structure, d'une veine de fluide 
polyphasique definie par plusieurs grandeurs physiques et comprenant des phases liquides 
et gazeuses. Suivant cette methode, on realise une modelisation du mode d'ecoulement : 

15 - en formant un reseau neuronal non lineaire avec une couche d'entree avec autant 

d' entrees que de parametres de structure et de grandeurs physiques necessaires a une bonne 
estimation de la sortie, une couche de sortie avec autant de sorties que de grandeurs 
necessaires a 1 'estimation du mode d'ecoulement et au moins une couche intermediate ; 

- en constituant une base d'apprentissage avec des tables predefinies reliant 
20 differentes valeurs obtenues pour les donnees de sortie aux valeurs correspondantes des 

donnees d'entree ; et 

- en determinant par iterations des facteurs de ponderation de la fonction 
d'activation permettant de relier correctement les valeurs dans les tables des donnees 
d'entree et de sortie. 

25 De preference, on analyse des donnees de sortie des neurones de facon a trier, parmi 

les valeurs des donnees de sortie du reseau de neurones, les seules donnees pertinentes a 
prendre en compte dans la determination iterative des coefficients de ponderation de la 
fonction d'activation. 
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La methode selon 1'invention 

- La methode selon V invention permet de construire un module (hydrodynamique 
ou thermodynamique parexemple) destine a simuler en temps reel le mode d'ecoulement 
en tout point d'une-conduite, d'une v.eine de fluide polyphasique comprenant au moins une 
5 phase liquide et au moins une • phase gazeuse, de fagon'qu'il. spit le mieux, adapts a des 
conditions operatoires fixees portant sur un certain nombre de.parametres structured et 
physiques definistelatifs a la conduite,.et-sur~un ensemble de grarideurs; physiques clefinies 
/^(grandeurs hydrodynamiques oil thermodynamiques par exemple), avec, des y gammes;de 
i < variation fixees pour les parametresiiet les grandeurs physiques: v >>/ 

/in >- Elle ''compoite rutilisWioil d'un systeme cle modelisatioh constitue a base v de \\ 

/ / c ;■ : g / F \ %\$. /S' \^J\ T • , V\. \ f > . \ ^ • 

// reseaux neuroriaux non lineaires ^e^ chacun^s fehtrees pour des parametres de structure V\ 
// ^" / k/ itiCT>^^^^^^v \\ \ \\ 

j ; et des grandeurs physiques, et des|s$Wretf 9U^ont / dispraibles des grand<Sjdrs\necessaires .a \ \ 

! f ?y ' \,kU(//'k , A . V\ \ T ,WW ; \\ 

f J ^estimation du mode d'ecoulemenf^vetXauomdins^unfe couGhej intermediaire^Les reseaux ■, 

il V ^^^^^M ,ul M 



l ! ' neuronauxr sont determines. :iterati^irient— points; aju&ter a 
I * 15 " d'apprentissage avetfdeOd}fe^ 



pour vs aiuster aux valeursi une base 



valeurs obtenues pour les 



I \ v -donnees helsbrtie^xivaeurs ^biTe^pondarites^es clonnees d' entree. / 

\\ ? 5ux\conditi0n\^operitolH^mpo5 optimises // 

\\ s'aiustantau mieuxkux conditions opfar6ires imposees. // / /'\ // 

X\ /V ^ // 

20 \\, ^Dans le cas, p"ar example, ou le module doit etre integre<a Un/modele g^neraha la Ms 
^fiydrodynanuque et thefmc^naS^ on 
utilise le mod^lelpour former la base^d-apprentissagerde maniere a selectipnner F ensemble 




: /ICS oils pmaiiiciicx\Ci: uiic^ giduu' 

25 Ton genere les reseaux^neuronaux optimises s'ajustant au mieja^a^a base d'apprentissage 



formee. ^^^- : -ZZ^~t:Z^^^'^' 
DESCRIPTION DETAILLEE 
Generalites sur la formation du systeme de reseaux 
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On considere une circulation de fluides polyphasiques dans une conduite avec au 
moins une phase liquide et au moins une phase gazeuse et Ton cherche a construire un 
sy steme neuronal permettant, a partir d'un certain nombre de donnees d'entree 
geometriques et physiques relatives a la conduite et de donnees physiques sur les fluides, 
5 de donner a chaque instant et pour chaque section de la veine fluide une estimation du 
mode d'ecoulement, ou bien, a partir d'un certain nombre de donnees liees aux conditions 
de pression et de temperature sous lesquelles se trouve le fluide, et a la composition du 
fluide, de donner le comportement thermodynamique du fluide. L'ensemble de ces 
donnees d' entree constitue une base de donnees. 

10 1) entrees et sorties 

Les donnees d' entree sont par exemple : 

- des donnees geometriques et physiques des conduites : diametre, rugosit6 et angle 
d'inclinaison du conduit, et des donnees thermodynamiques : masse volumique de gaz, 
masse volumique de liquide, viscosite du gaz, viscosite du liquide, etc. ; l'ensemble de ces 

15 donnees constituant un n-uplet ; et 

- des donnees hydrodynamiques caracterisant le melange : tension superficielle 
gaz/liquide, fraction volumique de gaz, vitesse barycentrique du melange, etc. 

ou bien : 

- des donnees thermodynamiques relatives a l'environnement du fluide : pression, 
20 temperature, l'ensemble de ces donnees constituant un n-uplet, ainsi que des donnees 

caracterisant la composition du melange : masse molaire, proprietes critiques, coefficients 
d'enthalpie, etc. 

Le modele produit par exemple en sortie le comportement hydrodynamique des 
effluents, et, notamment, le regime d'ecoulement. D evalue et d£livre sur deux sorties 
25 principals des donnees hydrodynamiques dans la partie de conduite dont on souhaite 
determiner le type d'ecoulement, la difference dV de vitesse entre gaz et liquide par 
exemple, ou la fraction (3 (j3e [0;l]) d'ecoulement de type stratifie (type ou la phase 
liquide s'ecoule dans la partie inferieure de la conduite). D'autres grandeurs qualifiant le 
type d'ecoulement peuvent etre calculees a partir de ces deux sorties. 
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Dans les applications a la thermodynamique, le modele produit en sortie les 
grandeurs thermodynamiques des effluents, comme, par exemple, le nombre et la nature 
des phases en presence ou les fractions molaires des phases. 

2) Structure du systeme de reseaux ; rf ^ , 

La fonction que doit model iser le modele hydrodynamiquef ou thermodynamique fait 
appela differentes sous-fonctions que l'on peutjesoudre en utilisant sLnecessaifeCplusieurs 
reseaux ddnt Te role au sein du [ .modele esTibieff defini.^Par exemple, !wv 'reseau peut\etre 

/ * v>*' . — rUiii^^ ^i .Js^ 




etc.). Ainsi, tin parlera par la suite de systeme; base suiydesrreseaux de neurones. 



* 1 



f ./d' entree 



15 r 



^referen< 
adapte a 
<fois les 



•encelde type-lNilEI^fe^ 



ri i ***** 

''art.tcaKil^est particulierement 'bien 

£ I// 

' 1 deerire toura-ra 



fixant meme^ si 



necessaire certams^vde^ses^parametres a. une^-valeurtf reifiee^efcfdonc physiquement 



w 

\\ representative A Ainsi, la connaissance d'equations physiques regissant les ecdQlements// 
2ti\\ permfet^par exemple ^d'enrichir le reseau et de l'adaptei^,au^mfeux aux <phenomenes v 
\p|iysique|^qijuil modehse is ^^ -^^^^^^ <Q 



// 



w 



// 



Le rese^uj/de neurones corripOrte~(-Figrl-)'^par exemple trois^,couches : la/couche 

-dessus) du 

spdndant aux deux 



25 



30 



d'entrdeCeomposee' de-^dix neurones correspondant aux dixfdonnees (citees/ci 
modele physique complet; une couche de sortie de^deux\neurones correspondai 
parametres recherch6s^l/et P, et une couche intermediaire^dite-'couche cachee, dont le 



nombre de neurones N c est optiMse^I^TreseaaSsWotalement connecte. La non linearite de 
ce reseau est obtenue par une fonction d' activation sigmoide regissant le comportement des 
neurones de la couche cachee. Les neurones de la couche de sortie peuvent etre choisis 
liheaires ou non-lineaires. La fonction d'activation peut etre par exemple la fonction 
identite. 
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3) Apprentissage 

Les poids de cette structure sont determines a Tissue d'une phase d'apprentissage ; 
au cours de cette phase, on nourrit le reseau d'un ensemble de donnees constituant la base 
d'apprentissage du reseau, et on optimise la structure et les poids du reseau en minimisant 
5 des erreurs constatees pour 1' ensemble des echantillons de la base, entre les donnees de 
sortie issues du calcul du reseau et les donnees attendues a la sortie, donnees par la base. 
Les erreurs peuvent etre les erreurs absolues entre les grandeurs d' entree et de sortie ou les 
erreurs relatives, selon la performance desiree pour le reseau. 

Les facultes de generalisation du reseau sont ensuite testees sur sa capacite a bien 
10 calculer les deux sorties pour des entrees qui lui sont inconnues. 

Particularites de mise en ceuvre 

Le systeme base sur des reseaux de neurones modelisant le module 
hydrodynamique ou le module thermodynamique est defini par un ensemble de parametres 
que la phase d'apprentissage permet de determiner et de fixer. Ces parametres sont 
15 adaptables aux contraintes fixees par l'utilisation a laquelle est dedie ce module. 

- Si le module est destine a une utilisation dans un cadre general, la base 
d'apprentissage utilise doit etre exhaustive c'est-a-dire inclure toutes les donnees 
disponibles concernant les differentes entrees/sorties ; le module est ainsi apte a representer 
Tensemble des cas que Ton peut rencontrer sur terrain. 

20 - Si le module est destine a etre utilise dans un cadre ou un contexte plus specialise, 

la production d'un champ petrolier par exemple, le calage est effectue sur une base 
d'apprentissage plus restreinte contenant seulement l'ensemble des donnees que Ton peut 
rencontrer sur ce champ uniquement : type ou structure des conduites ou pipelines, 
variations topographiques, gammes plus restreintes de variation de certaines grandeurs 

25 hydrodynamiques ou thermodynamiques, etc. Les parametres des reseaux neuronaux sont 
dans un tel cas calcules de maniere a obtenir un produit specifique a une configuration 
donnee, de fa<jon a representer au mieux les ecoulements dans des gammes operatoires 
propres a ce champ. 

Ainsi, ce module est adaptable aux differents objectifs fixant son utilisation. 



Environnement du module a reseaux neuronaux 

Le module developpe peut etre utilise seul, sans autres modules associes. Dans ce 
cas, il est cale sur des. bases de donnees constitutes en fonction de 1' utilisation desiree, 
selon que le cadre est general ou speeifique. II est independant et permet d'estimer les 
grandeurs hydrodynamiques ou thermodynamiques susceptiSles de permettre le suivi du 
comportement -des eqoulements d'effluents dans les pipelines, a ^partir tt'une base 
d-apprentissage regroupant des'ldonnies-de^^andeurs .physiques et^thenipdynamiques 
-fournies^par ailleurs. V . 
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Le module-pedFegalememietre irnpl*^ modele N (ensemble de modulesx 

logieiels de simulation capableside^imulei^le^cbmDortement transitoire\et/ou stationnaire \\ 

Am / Mr ^ ^ ^ % \ W 

des.ecoulements polyphasiques, tel ; lque^pap:*txeriiplg;;:4e modele TACfflE\deja cite^Dans \\ 

f $ mj}i%^f^. , v \ \ v - \\ 




\^ \ 

dpiinees d'emr^t^te^rtg^^ 

^^Cxemples d'util^sation 

i\\ lb, module a b^e de^reseaux neuronaux convientpour toute applicationnecessitant 
uhe representation hy<kodynarmq^^ de cpdes< de simulation 

% \\ ^ $ J) & — ^.•^-'^ ^ ^ y/ 

predictifs d'eeoulSments, de simulation temps reel ou de simulations' d'entrainement, etc. 

_ N ^ J „ _^si, ^ „ J. , ^L\r.._ yz, 



Les methodes de cotitrQle'actif-en production peuventaussbse^baser sur un.tel/module. 

X V v OF I a PW U " ^ 

Exemple>d!op^misation ^ &-/^ * ^ 

25 On considere par^xeinplerle^ un reseau donnant la 

classification du type d'ecoulement de fluides polyphasiques : ecoulement de type stratifie 
ou de classe CI (definie plus haut), ecoulement de type intermittent ou de classe C2 
(caracterise par une succession de bouchons liquides et gazeux), ou bien encore de type 
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disperse ou de classe C3 (le liquide etant entraine sous forme de fines gouttelettes), en 
fonction d'entrees comme decrit au paragraphe « Generalites sur la formation du reseau ». 

Dans un premier temps, on forme un reseau neuronal a l'aide d'une base generate. La 
structure optimale resultante que Ton obtient comporte 30 neurones, avec une seule couche 
5 cachee. La table de confusion du reseau obtenu (Fig.2) exprime le pourcentage de 
classification obtenu dans les trois classes de sortie du reseau pour chaque classe de la base 
de donnees. 

Dans un second temps, on restreint la base de donnees utilisee en apprentissage a 
quelques configurations particulieres, qui diminuent les plages de variation des differentes 
10 entrees du reseau. Le reseau neuronal optimal est dans ce cas constitue de 28 neurones et 
comporte une seule couche cachee. Le resultat de classification est materialist par la table 
de confusion de la Fig.3. 

On verifie immediatement que les resultats de bonne classification, repartis sur la 
diagonale, sont bien ameliores par rapport au cas precedent forme au moyen d'une base 
15 moins selective et moins adaptee a la configuration ou le module doit fonctionner. 

Ainsi, une adaptation de la base de donnees aux conditions d' utilisation permet de 
modifier d'une part la structure du reseau, d' autre part d'optimiser les resultats de 
generalisation obtenus. 
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REVENDICATIONS 



« 1) Methode pour former un module destine a__si muler en temps reel l e mode 
d'ecoulement en tout point d'u ne c oniduite , d'une vein6 de flui.de polyphasique comprenant 
5 au moins line phase liquide et au rrioins une phase gazeuse; de fa§on qu'il sqit le mieux 
adapte a des conditions operatoires .fixees portantsur un certain norpBre de parametres 
structurels ; definis relatifs a la^eoiiduite, et sur ''U^ejnseinbje de grandeurs physiques 
/.■ definiesfc avec des gammes de variation fixees' pouf les ^dits. 'parametres: "et les dites. 

/ ' ' - " j l v JPl ^%^}^A> " : ^ ; x\ N \\ 

// srandeurs physique#ku moyenJ^njsysteme^bH^^ consjitue & base de reseaux\\ 

f[Q neutonaux non lmeaires avec crtacuhl des eritr^^pp^des parametres ^ \\ 
// grifideurs physiques, et des sor&s^ des grandeurs \neces^aites y a \\ 

J restimation^u mode d'ecoulemeift|^i^^ les respaux \\ 

- neuronal^ |^tant deterpi^ valeu^ jfune bas£ 

; ^d'apprentisslage avec desvtabre? r ^ obteiiues pour lesj 

15 - donnees defsbrtie auxovaleurs ^ caractensee en ce que 

\ -Ton fornl |e#^|| 5 ^^ ! SQ£^|a»ires Jnjiosees 
\ ' . genere de^iiseaux neur(^^3pfifejes-s' ?ju.4^^M^^|^|^^o/s operatoires 

\\ Methods IJbn la revendication 1, caracterisee en ce queJ-'en^semble de grandeurs/ 

20\Xphysiqo^e,stconstitu^^^jdeurs hydrodynamiques. ^y^s // 

\b. 3) Methode selon la^enditl^ l'ensemblejde grandeurs 

\\ <' cr \ °' / / 

- - de prandeurs thermodvnamiaues. 

<0 



. , . .... „ r . , . aracterisee^en^ce que, ic un 

\n O F if A P K 1 ^ 7 

module etant integre4 un modele ^general a la fdis hydrodynamique^et thermodynamique 




25 de simulation d'ecoulements^p^phasj^ues, onj^tjlise^^tlirfnodele pour former la base 
d'apprentissage, de maniere a selectionner l'ensemble de grandeurs physiques le mieux 
adapte au fonctionnement du modele, ainsi que les gammes de variation fixees pour les dits 
parametres et les dites grandeurs physiques, et Ton genere les reseaux neuronaux optimises 
s'ajustant au mieux a la base d'apprentissage formee. 
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